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Abstract 



> 

in 
a^ 

' We investigate a relativistic self-interacting gas in the field of an exter- 

nal pp gravitational wave. Based on symmetry considerations we ask for 
those forces which are able to compensate the imprint of the gravitational 
' wave on the macroscopic 4-acceleration of the gaseous fluid. We establish 

I an exactly solvable toy model according to which the stationary states 

O ■ which characterize such a situation have negative entropy production and 

are accompanied by instabilities of the microscopic particle motion. These 
, features are similar to those which one encounters in phenomena of self- 

organization in many-particle systems. 

> 

X 



1 Introduction 

Let a gravitational wave (GW) act upon a simple relativistic gas which was 
at "global" equilibrium (see, e.g.,P El El before the infall of the GW. The 
influence of the non-stationary radiation field on the gas dynamics will generally 
destroy such an equilibrium. States of "global" equilibrium are only possible in 
stationary spacetimes, i.e., in spacetimes which admit a timelike Killing vector 
(KV) P [2^, _3j ^ ^ i6j . In the limiting case of ultrarelativistic matter (radiation) 
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the corresponding equilibrium conditions imply the existence of a conformal 
timelike KV. Neither a timelike KV nor a conformal timelike KV is compatible 
with the properties of gravitational radiation. A pp GW is connected with 
the symmetry groups G3, Gs and Ge, implying a covariantly constant null 
vector field and a set of spacelike KVs [71 The gravitational radiation will 
therefore induce non-equilibrium processes inside the gaseous system which can 
be described with the help of linear irreversible thermodynamics [51 110111X1112) . 

Typically, one studies the motion of particles or the behavior of continuous 
media under the influence of an infalling GW in order to reconstruct the prop- 
erties of the latter. This strategy relies on the circumstance that free particles, 
e.g., move on geodesies of the spacetime provided by the external field, and 
neighboring particles experience tidal forces which may be used to monitor the 
GW. Similar considerations hold for continuous media. 

The universality of the gravitational interaction makes the behavior of a 
material system under the influence of an external gravitational fleld different 
from the corresponding behavior under the action of an external electromag- 
netic field. In the latter case screening effects occur which may compensate the 
external action in the interior of the system. There is no corresponding shielding 
of the gravitational interaction. 

On this background one may, however, raise a different question: Is it 
possible to counterbalance (parts of) the imprint of an external GW by non- 
gravitational forces? Or, in other words, which kind of non-gravitational inter- 
action might (at least partially) "shield" the action of the gravitational field? 

In the present paper we investigate this problem by admitting suitable self- 
interactions of a Maxwell-Boltzmann gas as candidates for such hypothetical 
screening. The latter will be characterized by symmetry considerations as fol- 
lows. Before the GW infall the system is assumed to have certain spacetime 
symmetries, e.g. those of Minkowski space, which, in general, are broken by the 
GW. Our requirement for a (partial) compensation of the GW imprint on the 
gas will be that a part of the original symmetries continues to hold even under 
the GW influence. Then we show, how this requirement can be realized with 
the help of the mentioned interactions. 

This strategy is comparable to the one used in gauge theories. Gauge field 
theories rely on the fact that local symmetry requirements neccessarily imply 
the existence of additional interaction fields, the gauge fields. Although the 
symmetries in this paper are of a different type, the general feature, that sym- 
metry requirements are accompanied by additional interactions, is true here as 
well. 

On the other hand, there is no guarantee, that the required forces really do 
exist. In case they do not exist, there is no shielding of the described manner. 
But in any case, establishing the self-interaction which would be necessary to 
realize a given symmetry may improve our understanding of the dynamics of 
physical systems in external gravitational fields. 

This paper relies on a description of the gas by the relativistic Boltzmann 
equation for the invariant one-particle distribution function. According to the 
conventional picture the gas particles follow geodesic worldlincs of the radiative 
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spacetime in between elastic binary collisions. Imposing symmetry conditions 
of the mentioned kind will, in general, make the particle trajectories differ from 
geodesies. This difference to geodesic particle motion may be attributed to 
interactions which can be described either by an additional (in general inelastic) 
collision term or by an effective force on the individual particles. 

The crucial point is that this interaction is not known in advance. It depends 
on the external field and represents that response which is necessary to realize 
the imposed symmetry. It is equivalent to a self-interaction of the gas. Either 
there exist corresponding interactions inside the gas or these forces are mimicked 
by an external source. Both options are treated here on an equal footing. If the 
system is unable to respond in the required manner, the symmetry cannot be 
maintained. In other words, the system has to organize itself properly in order 
to counterbalance the external field successfully. 

The considerations of this paper represent an exactly solvable toy model 
along these lines. We determine the necessary interaction for a special case. As 
an initial configuration we take a homogeneous and isotropic gaseous fluid in 
Minkowski space which macroscopically is at rest. Under these conditions it has 
a vanishing 4-acceleration — 0. In the field of a subsequently infalling GW 
the fluid motion will no longer be geodesic in general, i.e., one has f/* 0. In the 
present case the mentioned symmetry requirement will be equivalent to enforce 
a macroscopic geodesic fluid motion — also under the influence of the GW 
field. This is achieved by allowing the trajectories of the individual microscopic 
gas particles to be non-geodesic. (Without imposing the additional symmetry 
it would be vice versa, i.e., the microscopic trajectories were geodesic, but not 
the macroscopic fluid motion.) The corresponding states are stationary, far- 
from-equilibrium states, accompanied by instabilities in the microscopic particle 
motion. There exist striking similarities to processes of self-organization |13II14[ 

nsj. 

The paper is organized as follows. In Section |2] we determine the effective 
one-particle force on microscopic gas particles which is compatible with an equi- 
librium distribution of the latter (generalized equilibrium). We establish our toy 
model by suitable symmetry considerations concerning the vector U°' /T, where 
If^ is the macroscopic 4- velocity and T is the temperature. In Section O we 
derive the corresponding macroscopic balance equations for the first and sec- 
ond moments of the one-particle distribution function which, as a consequence 
of the force action, have nonvanishing sources in general, characterizing the 
self-interaction of the system. In Section 0] the equations of motion of the mi- 
croscopic particles are explicitly integrated. The self-interaction is shown to be 
accompanied by instabilities of the transverse particle motion. Our main results 
are summarized in Sec. |S1 Units have been chosen so that c — ks = h = 1. 
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2 Force field and generalized equilibrium 



The one-particle distribution function / = f {x,p) of a relativistic gas is sup- 
posed to obey the Bohzmann type equation 

L if] EE p^f,, -Tlip'^pi^^ ^J[fJ]+S ix,p) , (1) 

where J[f, f] is Bohzmann's colhsion integral. Anticipating that additional in- 
teractions will play a role in the following, we have included here a non-standard 
term S on the right-hand side, which takes into account those "collisions" which 
can not be reduced to elastic, binary interactions. Following earlier applications 
m a cosmological context [TBI [T71 [THIIT^ we assume S to have the structure 

It is straightforward to realize that a "collision" term of this form may be taken 
to the left-hand side of Boltzmann's equation resulting in 

P'fn-niP'p'§,+mF^§^^-J[fJ] . (3) 

The left-hand side of this equation can be written as 
df{x,p) _ df dx' , dfdp' 



d'y dx^ d'y dp^ d'y 

with 

d7=^' dT^'"^- 

Equations I^J are the equations of motion for gas particles which move under the 
influence of a force field — {x,p). The quantity 7 is a parameter along the 
particle worldline which for massive particles may be related to the proper time 
T by 7 = r/m. Consequently, a specific "collisional" interaction, described by a 
"source" term S, may be mapped onto an effective one-particle force F'^. This 
demonstrates that there exists a certain freedom to interpret collisional events 
in terms of forces. (This freedom can also been used in the reverse direction, 
i.e., to interprete (parts of forces) as collisions "W). We emphasize that our 
approach is different from the "canonical" theory of particles in a force field for 
which the force term mF^df /dp^ in Eq. ^ is replaced by md {F^f) /dp^ |2()| . 
While both approaches are consistent with the equations of motion Q, they 
coincide only for dF'^/dp^ — 0, which holds, e.g., for the Lorentz force. In the 
cases of interest here we will have dF'^/dp'^ ^ 0. 

For collisional equilibrium / reduces to Jrittner's distribution function Z*^, 

/^/°(2;,p)-exp[a + /3,pT , (5) 
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where a — a(x) and /3a ix) is timelike. For / ^ /° the colhsion integral 
vanishes, i.e., J [/", /"] = 0. Substituting f into Eq. Q we obtain 



(6) 



As already mentioned in the introduction, the situation here is different than for 
example for charged particles in an electromagnetic field where the microscopic 
starting point for a self-consistent analysis is a given interaction, the Lorentz 
force. Here, the analytic structure of the force is not known in advance but has to 
be determined. This requires a kind of "inverse" strategy which is characterized 
by symmetry considerations. In the present paper we are interested in that 
specific force which realizes the requirement V^ffsgik = where /3* = /T. 
This corresponds to retaining a part of the original timelike symmetry fpgik = 
0, by which the system was characterized before the GW infall. 

Once the analytic structure of the force has been determined on this basis, 
the analysis may be carried out along the usual lines of a self-consistent treat- 
ment. Focusing on forces which satisfy the equilibrium condition ®, the most 
general structure of the term on the right-hand side of equation © is 



(7) 



Comparing now different powers in Eq. © separately, we obtain the (general- 
ized) equilibrium conditions 



and 



kl 



(8) 



(9) 



The last equation relates space-time symmetries to the force exerted on the 
particle. The Lie derivative may be invariantly split into contributions parallel 
and perpendicular to the fluid four- velocity: 



where hi 



9ik 



£p9ik — —UiU"' £ pgnk + h^£p9nk , 
U,Uk, 



(10) 



and 



U^£0gnk - --^^Uk + ^ 



hi£09nk — — 



Uk 
T 



T 



Uk 



2(7. 



^kT 



T 



ik 



T 



29 

1 

3T 



ik • 



(11) 



(12) 



Here, T = T^aU"" is the derivative in direction of the four-velocity, Ua = Ua-cU'^ is 
the four-acceleration, VaT = h^T^t, is the covariantly defined spatial derivative, 
Q = U\ is the fluid expansion scalar, and aab is the shear tensor. The quantity 
V/T is a Kilhng-vector (KV) if both parts (HH) and in (HHJ are zero. 



5 



Let us assume that the initially homogeneous and isotropic gas is at rest and 
in equilibrium and that £pgik — holds before the GW infall. Generally, this 
symmetry will be destroyed by the GW. The action of a GW on matter is known 
to produce anisotropies and inhomogeneities within a test medium of this kind. 

Now we ask, whether and under which conditions it is possible to retain a 
part of the latter symmetry, namely the projection 

u'fpgik = , (13) 

also in the GW field. According to this condition implies T — Q and Um + 
^^"^ = 0. This is equivalent to retaining that part of the previous equilibrium 
conditions which fixes the temperature behavior but to allow the quantities G 
and (Jab to be different from zero. This case turns out to be exactly solvable as 
will be shown in the remainder of the paper. For the special case VmT = the 
above condition implies Um = 0, i.e., a geodesic fluid motion. 
Imposing the corresponding "partial symmetry" condition 

2 9 2 

£ fjgik ^ Q £ (jgik = -^Tj^^ik + —Uik 

1 2Qik 
= ^[V,C/, + V,C/fe] = ^ , (14) 

is consistent with a non- vanishing quantity UiF^ in which 

mU.Fl^ = , (15) 

i.e., this condition fixes the part which is quadratic in the particle momenta. 
Obviously, macroscopic fluid quantities enter the analytic expression for the 
force which is exerted on the microscopic particles. The effective interaction 
couples the microscopic particle momenta with quantities such as the fluid shear 
tensor which characterizes the system on a macroscopic level. This is quite 
similar to the well known self-consistent coupling of the particle momentum of 
a charged particle to the electromagnetic field strength tensor. The linear part 
of the force is determined by the derivative of the quantity a according to 

Tfl 

- yC^.^a = -aC^a+Vaa . (16) 

The partial symmetry condition H14|l and Eq. H16|l constitute the conditions of 
"generalized equilibrium" for the system of gas particles in the force field Q . In 
the force-free limit these conditions coincide with those for "global" equilibrium 

PEHH. 

An example for a force linear in is the Lorentz force which is obtained 
for F"^^ m-^iem) ^hcrc e is the charge and F^^^^^ is the electromagnetic field 
strength tensor. For this case the equilibrium condition © is well known in 
the literature |21 E] . Other cases have been discussed in a cosmological context 

nsnnmanni 
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The conditions (|15|l and do not yet fix the spatial part of the force. The 
condition for the latter is that the force itself has to be orthogonal to the four 
momentum, i.e., PiF^ = . This requirement is satisfied for the choice 

mF^.pV = [QabU' - eiUb] . (17) 

For particles moving exactly with the macroscopic fluid velocity, i.e., p^ — mU^ , 
this force vanishes. It follows that the force is related to the relative motion of 
the individual particle and the fluid as a whole. In a sense this is similar to the 
motion of a sphere in a viscous fluid described by Stokes' formula. The spatial 
part of this contribution to the force is 

mK:Fl,py^ = e;'i?p" , (18) 

the quantity E = —UaP'^ being the energy of the particle for a comoving ob- 
server. In the nonrelativistic limit with E fa m we obtain 

h^^Fl^p'^p" ^ 1 (i;,,^ + v^.,)p^^ , (19) 

where is the nonrelativistic fluid velocity and /x, i/ = 1, 2, 3. This force couples 
the particle velocity p^ jm to the velocity derivative of the medium through 
which the particle moves. It represents a kind of friction. In this analogy the 
symmetrized spatial derivative of the fluid velocity plays the role of a friction 
"constant" . The difference to the motion of a nonrelativistic particle in a viscous 
fluid is the self-consistent nature of the problem in our case: The particle which 
is subject to the force exerted by the medium is part of this medium itself. 
In the limit of Eq. H19|l this force represents a linear coupling of the particle 
momentum and the derivatives of the fluid velocity. In a sense, this parallels the 
Lorentz force which couples the particle momentum linearly to the derivatives 
of the electromagnetic potentials. 

Finishing this section we emphasize again the unconventional strategy which 
was used here to obtain the analytical properties of the force F"^: The latter 
quantity which is similar to a friction force characterizes just those interactions 
which guarantee the prescribed (partial) symmetry condition (|14|l of our model. 



3 Transport equations 

The particle number flow 4-vector A^* and the energy momentum tensor are 
defined in a standard way (see, e.g., 013) as 

7V^ = y dPp7(a;,p) , r^^^y APfp^!{x,p) . (20) 

The integrals in the definitions H2U|) and in the following are integrals over the 
entire mass shell p^pi — ~rr? . The entropy flow vector 5° is given by [2121 

S- = - lp-[flnf-f]dP , (21) 
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where we have restricted ourselves to the case of classical Maxwell-Boltzmann 
particles. The entropy production density is a sum of two terms: 

Si^^(7c + (JF . (22) 

Here, 

<J, = - J dPJ[f,f]\nf (23) 
is the familiar contribution of Boltzmann's collision integral, while 

ap=m j APF'^Anf (24) 

takes into account an entropy production due to the action of the force F"^ . 
Since Boltzmann's H theorem guarantees Uc > 0, we have 

S^i-aF>0 . (25) 

The equality sign in the last relation is realized for f ■ With / replaced 

by /° in the definitions and the quantities 7V°, T"^ and S°- may be 
split with respect to the unique 4-velocity C/" according to 

N^^nU", T'^ ^ pV'U^ + ph"^ . S^^nsU''. (26) 

The balances 



may be written as 
and 



T-;^ = -mA / p"-F'fdP , (27) 



TV.: = -mf3^ {F^N^ + F^.T^") , (28) 

T"^ = -mP, {F^T'"'' + FliM"^^) , (29) 

where M = / AP f^p°^p*'p^ is the third moment of the equilibrium distribution 
function. For the entropy production density we find 

S'^, = mA y [a + /3aP°] F7"dP - -aTV.^ - /J.T'^g . (30) 

Evaluating the integrals in the balances (|28|l and (|29|) we obtain 

iV° ^na + nQ = nT , (31) 

UaT''.!l^-ap-{p + p)e , (32) 

and 

hcaT'^l = pVca . (33) 
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Before the GW infall the gas was supposed to be at rest and the initial timeUke 
symmetry imphed n = T — p — s — 0. Now we ask under which conditions the 
quantities n, T, p and s may be constant in time also under the GW influence. 
Obviously, this requires a force which guarantees d = 0. Indeed, as to be seen 
from Eqs. and (|32|l with (|26|l . the latter condition implies n = p = 0, while 
T = is already part of our basic assumptions [cf. the discussion preceding Eq. 
H14() ]. From the Gibbs equation 

Tds = d- + pd- ^ nTs = p-[p + p)- (34) 
n n n 

it follows, that s = is valid as well. However, the force which maintains the 
equilibrium distribution gives rise to source terms in the balances H31|l - (I33|) . 
Two interpretations are possible here. Either the sources (right-hand sides) 
in H32|l and (|33|l may be taken to the left-hand sides of these equations and 
regarded as stress components of a conserved energy-momentum tensor for a 
self-interacting medium [cf. 1171 IT51 or they describe an equivalent, 
but externally driven interaction which effectively is a self-interaction as well. 
Dynamically, both interpretations are indistinguishable. The interaction terms 
describe the support from internal or external sources which is necessary to resist 
the "attack" of the GW. However, our toy model leaves open the question, how 
these interactions my be realized. For the part ap of the entropy production 
[cf. Eq. (|24|) ] we find under these circumstances 

(jp = nsQ , (35) 

where we have used the expression 

P+P 



nT 



(36) 



for the entropy per particle together with the identification a — p/T where 
/Lt is the chemical potential. We emphasize that our generalized equilibrium 
due to the partial symmetry requirement (|14|l is consistent with a nonvanishing 
entropy production density. While the contribution cjc [cf. Eq. H23|l ] to the total 
entropy production density (|22|l vanishes for f f'^, the contribution ap [cf. 
Eqs. H24|l and l|35() ] does not. The corresponding states are out-of-equilibrium 
configurations of a specific kind. Generally, the theoretical description of out-of- 
equilibrium phenomena is known to be a major problem as soon as the deviations 
from a standard equilibrium reference state cannot be considered as small. The 
concept of generalized equilibrium is an approach to characterize a certain class 
of states which may be far away from "conventional" equilibrium (vanishing 
entropy production), but nevertheless are described by a Jiittner function and, 
consequently, admit an exact analytic treatment. Macroscopically, the imposed 
(partial) symmetry implies 



S^i, = nse , h = p^s^O . (37) 

This means that the equilibrium properties ((Sj and n = T = p = 3 = 0, 
which are supposed to characterize the system before the GW infall, remain 
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valid also under the influence of the gravitational radiation. It is obvious from 
Eq. I31|l that for such a situation to be realized, a change in the number of 
particles is necessary. The existence of source terms in the balances for the 
moments of the distribution function suggests an interpretation of the gas as 
an "open" thermodynamical system. It is well known that stationary states 
of open systems away from "conventional" equilibrium are key ingredients for 
phenomena of self-organization (see, e.g., ^J^^E] and references therein). 



4 Particle motion 

4.1 Energy conservation and expansion scalar 

In this section we study the microscopic particle motion in the force field H17|l 
under the influence of an external gravitational wave. According to Eqs. 
0, ((T7|) and neglecting the linear part of the force, i.e., PiF^ = we have to 
solve the system 

[QkiU, - QuU,]p^p' ,^ = p!: (38) 
L>T m clT m 

in the background metric 

dr^ = -2dudv + gaf3(u)dx°'dx^ , (39) 

where ^ ^ 

a,/3 = 2,3, u^—j=-, v=—j=-. (40) 

This metric is an exact solution of Einstein's vacuum equations with a G5 sym- 
metry group (see, e.g., [7], §33, |H1, chapter 21.5). It admits flve Kilhng vector, 
three of which, namely, 

CW='^- C[2)='^2, ei3)='^3, (41) 

being the generators of an Abelian subgroup G3. The flrst of the vectors H41|l 
is covariantly constant and isotropic, ^^2) '^(3) spacelike. Evidently, the 
dynamic equations H38I) admit the quadratic integral of motion = — m^. To 
flnd another integral of motion let us realize that 

^(t/>,) = v'-Vr,i + V'F, = -lpVc/,;7, . (42) 
UT m m 

In order to make further analytical progress, we first realize that the KVs 141|l 
of the radiative spacetime are also part of the symmetries of the initial config- 
uration before the GW infall. Under this condition it is tempting to assume 
that the fluid quantities "inherit" the spacetime symmetries |12l I21j . In partic- 
ular, we require that the symmetries £^^^^gik — for r = w,2,3, corresponding 
to the Killing vectors (|41|l . are inherited by the temperature, i.e., £^^,.-^T = 0. 
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Combining this kind of temperature behavior along the vectors (I41|) with the 
equilibrium condition T = implies a constant temperature both in space and 
time. According to the condition [/„ + ^-fr- = 0, however, a spatially constant 
temperature is equivalent to Ui = 0. Consequently, a combination of the con- 
dition U''£f3gik = with the requirement £^^^^T = of symmetry inheritance 
implies a geodesic fluid motion. Under these circumstances we have 

T>E 

— = , E = -U^p" = const , (43) 
ut 

i.e., the quantity E, which may be interpreted as particle energy for a comoving 
(with J7*) observer, is a second integral of motion. Contracting f/j with the 
Killing vectors H41() leads the following well-known consequences [2]: 

= {e(r)Ui)-=Q, (44) 

equivalent to 

e^,)f/, = const , Ur = Ur{Q) . (45) 
It follows that the unit time-like velocity vector in the GW field has the form 

C/, = S'iUM + SfU2{0) + 6fU3{0) + 5lUu{u) , (46) 

where 

^ . (47) 

This solution, which satisfies the symmetry inheritance condition £^^^^ Ui — 
for r = w, 2, 3, allows us to calculate the expansion scalar 9: 

e = ^^^JU-V^)=-U.m^JnV^,. (48) 

Here, = ~Uy is a non- negative quantity by definition and —g, where g is 
the determinant of the metric tensor, is monotonically decreasing with u (see, 
I22j §35.9). We find that 6 < 0, i.e., the expansion scalar is negative, 
which makes the present situation completely different from previous studies of 
generalized equilibrium in a cosmological context ^1 117L 1181 [TO] . According to 
the balance equation H31|l (with d = 0), a negative expansion scalar implies the 
annihilation of particles. As a consequence, we have negative entropy production 
as well [cf. Eq. H37|l ]. We emphasize that this kind of behavior does not 
contradict the second law of thermodynamics since the gas under consideration 
here is an open system (see the discussion following Eq. (|S3)- Loss of entropy 
in an open system is typical for processes of self-organization '13' . We argue 
that the mobilization of suitable forces to maintain the relations lO and 1)3 7|l 
in the presence of the GW exhibits features of a self-organizing system. The 
situation may be illustrated by a martial picture: Under the "siege" of the GW 
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the system will "evacuate" part of its population in order to maintain partially 
its previous internal state. Unless the GW stops to act upon the system this 
process continues until the system has lost all its particles. On the other hand, 
if the system is unable to organize suitable supporting forces, it will not succeed 
in further following its own internal dynamics but it will be determined by the 
external field. If there are no forces mobilized at all, the system is entirely 
"slaved" and its particles have to move entirely as prescribed by the external 
agent, i.e., along the geodesies of the radiative spacetime. 



4.2 Longitudinal motion 



By using relations (|43|) and 147|) the components p„ = —p" and py 
be expressed in terms of Pa and the components of Uf. 



Pv 



1 



e± 



y^£2 _ (1 + g'-^U^Uff) (m2 + g'-^p^pp) 



-p may 



(49) 



with the abbreviation 



and relation 



e = E~ g'^^Uo^pp 



2puPv 



g^^PaPlS 



(50) 



With the well-known interpretation of the covariantly constant isotropic Killing 
vector ^^^^ as wave vector of the gravitational wave (see, e.g., |23j, chapter 9, 
1211, chapter 15.3), the projection p„ = _Pi^(^-, of the particle momentum onto the 
GW wave vector represents the "longitudinal" component with respect to GW 
null direction. This scalar quantity will play an important role in the further 
integration process. The signs minus or plus in the expression H49|l correspond 
to the cases where the particle momentum has positive or negative projection, 
respectively, onto the spatial propagation direction of the gravitational wave. 
This may be most simply demonstrated if the system as a whole is at rest, i.e.. 



U2 = U^ = 0, Ui = 



Uv~Uu 

V2 



= , Uu = = - 



1 

V2 



(51) 



and consists of massless particles, moving parallel or anti-parallel to the GW 
propagation direction (z^-direction) with pa — 0. Under these circumstances 
we obtain 

(52) 



Pv 



-p 



-p"+p^ -EtE 



^/2 ^/2 V2 

The projection py vanishes for a positive p^, which corresponds to the minus 
sign in formula (|49|l . The plus sign in formula (|49|l together with the fact that 
Uu is negative according to Eq. H51|) . corresponds to a negative projection py, 
consistent with Eq. (151: 
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4.3 Transversal motion 



According to Eqs. H38|) and H43|l for Ua = 0, the components Pa satisfy the 
nonlinear equations 

dpa E 1^ 

-i- = —OakP ■ (53) 
dr m 

Their evolution is coupled to p^, and p„, which are given by formulas (|49|l and 
H5U|I . respectively. In order to simplify the analysis, we assume that the gas was 
macroscopically at rest before the infall of the GW. The propagation direction 
of the GW will be the only preferred direction in such a case and, we will refer 
to the components of the particle momentum perpendicular to this direction as 
"transversal" . If so, the relevant components of the tensor Qik are 

Qau ^eav=0 , QaP = ^^g^^(w) , (54) 

where the prime denotes the derivative with respect to u. Substituting the 
derivatives with respect to affine parameter r by 

_d_ _ du_d_ _ p^_d_ _ _d_ 
dr dr du m du to dw ' 

we obtain from Eqs. 149|l . 1)50(1 . and 1)53(1 the following two-dimensional nonlin- 
ear, nonautonomous dynamic system: 



Ep^ 



E± y/E^ -m^ - g"f^{u)pap(3 
A special solution of the last equation is 



(56) 



Pv = 



E ± y/E'^ - m? E T ^/E'^ - m 



2 



V2 ' V2 



p„ = , pi = t/b2~^ . (57) 

All components of the momentum are constant in this case, which corresponds 
to a vanishing of the self-interacting force on the right-hand side of Eq. 1(53(1 
with the tensor components 1(54(1. The solution ((57(1 is consistent with Eqs. ((49(1 . 
l|Sn|) and 



4.3.1 Linear stability analysis 

In a next step we consider the stability of the solution l(57() in the vicinity of 
Pa = 0. To this purpose we linearize the dynamic system 1(56(1 with respect to 
small deviations from zero {\pa\ << m). For simplicity we assume that 1723 = 0. 
Neglecting quadratic terms in pa we find, that the expressions 1(57(1 for p„ and Pv 
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remain valid up to first order in the deviations from pa = 0. The corresponding 
solutions are 

P2{u) = C2|.922Hr(^^ , PaH = CalssaHr^^' , (58) 

with 

, ^ E(Et VE^ - m2) ^ ^ 

nE) ^ > . (59) 

Obviously, r{E) is positive definite. Recalling that the lower sign in the last 
expression corresponds to a particle motion parallel to the spatial propagation 
direction of the GW, T{E) increases monotonically from i for = m, to (■^) , 
for E >> m. If the particle moves in the opposite direction, r{E) decreases 
monotonically from i for £' = m, to j for E >> m. The stability properties 
depend on the values of |522(u)| and |g33(u)|. In the absence of the GW one has 
|ff22(M)| = \g33{u)\ = 1 withp2(w) = C'2 andp3(M) = C3, independently of r(E). 
In the GW field, one of the functions [322 (m)| or 1333 (m)| is larger, the other one 
smaller than 1. For example, let us consider the well-known exact solution of 
Einstein's equations, where 1722 — cos^ ku and 1733 — cosh^ ku. (This solution is 
obtained from formula (25.23) in [7] by substituting c = c — 0, a — —1, = ku 
and = v/k, and by changing the metric signature.) The perturbations of the 
component P2 will decrease in this case, while there is a dynamical instability in 
the P3 component. The corresponding growth rate depends on the energy E and 
on the direction of the particle motion with respect to the spatial propagation 
direction of the GW. We conclude that the self-interaction which is necessary to 
maintain the (generalized) equilibrium under the influence of the external GW, 
is connected with an instability of the transverse particle motion. Dynamical 
instabilities are known to accompany processes of self-organization (see, e.g., 
[I'M 1141 ITS] '). Although a more detailed correspondence to such kind of effects 
in synergetics has still to be worked out, the existence of instabilities in the 
particle motion seems to support again the interpretation that self-interaction 
in the present context exhibits features of self-organization. 

4.3.2 Exact one-dimensional solution 

With the above result of the linear perturbation analysis that the particle motion 
in direction is unstable, while there is stability along the x'^ axis, it seems 
reasonable to investigate a model with p2 = 0, but p3 7^ 0. The equation to be 
solved in this case is 



dp3 E 533 ' 
du 2 333 



E± 



{P3 



933 



(60) 



\P3\ ^ cosh ku J E^ - - (Z - E)\ (61) 



Using again the solution (733 — cosh^ ku and substituting 
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the function Z is implicitly given by the solution of the transcendent equation 



[A-Z + ii;]«+MA + Z-£;]-*+^ = ^^^ , (62) 

cosh ku 

where denotes the abbreviation = ~ m? . An exact analysis is possible 
for massless (to — 0) or ultrarelativistic (to 0) particles. In this case one 
obtains 

_ „ const 

Z = 2E — , 63 

cosh ku 

and 



pi = Ci,+ 2E{coshku -1)1 E± JE^ - C| . (64) 



The integration constant has been chosen such that ps — C3 a.t u — 0. The 
trivial solution p3 = of Eq. Ht)U|) is realized in the expression (|64|) for C3 = 
with the minus sign in front of the square root on the right-hand side of 
this equation. The exact, explicit solution H64() shows that ps really increases 
with the retarded time u. It confirms the result of the previous linear analysis 

according to which the transverse particle motion is unstable. It should be 

2 

mentioned, however, that the quantity ^^^^^ , decreases, so that the expression 
under the square root on the right-hand-side of formula H6()|l remains positive. 



5 Summary 

Given that a homogeneous, isotropic fluid at rest in Minkowski space admits 
a timelike Killing vector [/"/T, we studied the behavior of this system under 
the condition that the projection £pgik of the Lie derivative of the metric 
in direction of remains zero also in the GW field, i.e., U'^£pgik — 0. This 
restricted timelike symmetry requires the existence of interactions which we 
have modelled as a self-interacting effective one-particle forces that is exerted 
on the microscopic particles of a Maxwell-Boltzmann gas. As a toy model 
we investigated the exactly integrable special case in which these forces result 
in compensating the GW imprint on the macroscopic fluid acceleration. We 
suggested an interpretation of the gas as an open thermodynamical system. In 
particular, the considered interaction is accompanied by a reduction of the phase 
space which is equivalent to a decrease of the entropy content of the system. 
The corresponding states are stationary and far from equilibrium. Moreover, 
they are characterized by dynamic instabilities in the transverse (with respect to 
the propagation direction of the GW) microscopic particle motion. We pointed 
out the similarity of these properties to the phenomenon of self-organization. 
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